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Abstract 

The B c ( 1 5o) meson to S-wave Charmonia transition form factors in large recoil region are cal- 
culated in next-to- leading order (NLO) accuracy of Quantum Chromodynamics(QCD). Our results 
indicate that the higher order corrections to these form factors are remarkable, and hence are 
important to the phenomenological study of the corresponding processes. For the convenience 
of comparison and use, the relevant expressions in asymptotic form in the limit of m c —> are 
presented. 

PACS number(s): 12.38.Bx, 12. 38. St, 12.39.Hg, 14.40.Nd. 
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I. INTRODUCTION 



The study of B c meson is of special interest, since it is the only heavy meson composed 
of two heavy quarks with different flavors. The B c exclusive decays provide an important 
non-relativistic system in the investigation of weak interaction, hadronic properties of heavy 
mesons and even new physics. Up to now there are only two decay modes of B c meson being 
observed in experiment at the Fermilab Tevatron, i.e. B^Sq) — > J/ipn and 5 c ( 1 5'o) — > 



J/ipe + v e [1]. Theoretically, many of works were carried out in different frameworks, see for 



instance recent works 



0,Q 



and references therein. In analyzing the B c decay processes, there 
are several different scales should be taken into account: the hard scale set by the heavy 
quark masses trig, the soft scale set by mqv where v < 1 is the relative velocity of heavy 
quarks within the B c meson, and the ultrasoft scale set by mqv 2 . The hard part supplies 
the short- distance contribution and can be calculated perturbatively in strong interaction, 
while the soft and ultrasoft parts belong to the long-distance contribution and have to be 
evaluated via some non-perturbative methods or fitted by experimental data. 



In the study of i?-meson decays, the factorization (4, 5] is crucial to disentangle the 
short-distance sector from the long-distance sector, where the former can be treated by 
perturbertive QCD(pQCD), while the later can be characterized by some universal hadronic 
parameters. Because B to light hadron exclusive decays are mediated by weak interaction, 
it is convenient to use an effective weak Hamiltonian to describe the interaction, which has 
the following structure: 

^ff = %Y, V CKMCi{^)Qi. (1) 
* i 

Tere Gp is the Fermi constant and Qi are local operators, Ci are short-distance coefficients 
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and Vqkm * s CKM matrix element 



9] . In naive factorization approach, the B- meson 



exclusive two-body decays can be formulated as 

{M 1 M 2 \Q i \B) ~ (M 1 \^Tb\B){M 2 \^TiP\0) , (2) 

where the matrix element (Mi\ifjTb\B) stands for the transition form factor at large recoil, 
and (M2\ipTifj\0) corresponds to the M 2 decay constant. If we consider all the partons 
on light-cone, the matrix element {Mi\ipTb\B) can not be factorized further due to the 
divergence at the end point arising from vanishing energy of the partons on the light-cone. 



However, if we extend to the non-relativistic situation, the nonperturbative effect can be 
factorized to Coulomb potentials of initial or final bound states. For instance, in the process 
B c — > J/ip(r] c ) + 7T, we can describe the dynamics of bound states B c and J/tp(i] c ) by non- 
relativistic QCD(NRQCD) 10], since the masses of bottom and charm quarks are much 



bigger than Aqcd- And then, the matrix element relevant to the form factor at large recoil 
can be factorized as [ll, 12^ : 



(3) 



Here, the nonperturbative parameters ipB c {fy and V'JMfac)^) are the Schrodinger wave func- 
tions at the origin for be and cc systems, respectively. T, is a hard scattering kernel which 
can be calculated perturbatively. 

As the LHC will soon be in the position to explore many B c decay channels - among these 
several are in semileptonic and nonleptonic charmonium decay modes - a dedicated study 
of 5 c -to-charmonium form factors is meaningful. In this work, we will explicitly calculate 
the matrix elements (J/ip\cTv<A)b\B c ) and {ri c \cTyb\B c ) in pQCD approach at the next-to- 
leading order in non-relativistic limit of the initial and final bound states. The paper is 
organized as follows: in section II, we represent matrix element at the Born level; in section 
III, we calculate the matrix elements in the NLO accuracy in pQCD; in section IV, we 
compare result from pQCD with the wave-function overlap contribution qualitatively; in 
the last section a brief summary and conclusions are given. 

II. THE FORM FACTORS AT BORN LEVEL 
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FIG. 1: The leading order Feynman diagrams 



The investigation of process B c decays to S-wave charmonia {J/ip or i] c ) with a light 
meson or lepton pair plays an important role in the study of B c property, where the nature 



of the transition form factor stands as a central issue. In this work, we focus on the study 
of two and four independent form factors in B^Sq) to i] c and J/ip transitions respectively, 
which are normally defined as: 

( Vc (p')\^b\B c (p)) = f + (p' + py + f_(p - p'Y , (4) 

(J/^(P',e*)\cYb\B c (P)) = ige^elP a P' p , (5) 
(J/ip(P', e*) \c^ 5 b\B c (P)) = a e^ + a+e* • PP'» + a_e* • PP^ . (6) 

At the leading order in a s , there are two independent Feynman Diagrams for 
(J/ip\cT v ^A)b\B c ) and (r) c \cr V ( A )b\B c ) , as schemetically shown in Figure 1. In non-relativistic 
limit, the momenta of constituent bottom and charm quarks are Pb = £,P and p 5 = (1 — £)P 
with £ = m ™£ for B c meson, and p 5 = p c = P'/2 for J/ip or i] c meson. Here, P and P' 
signify the momenta of initial B c and final charmonia. 

After taking the above mentioned procedures, it is straightforward to calculate those 
concerned form factors at the tree level. They read 



j, LO _ 8V2TTa s i/j(0)B c ip{0) Vc C A C F (^m b + m c )(3ml + 2m c m b + 3m 2 c - q 2 ) 

N c mJ' 1 '(ml + m 2 — 2m b m c — q 2 ^ 2 



j-lo _ ^QV^7ra s ip{0) B c ^{0) Vc C A C F {m b + m c ) 3/2 (m b - m c ) 
NjnJ 2 (ml + m 2 — 2m c m b — q 2 ^ 2 



LO _ 32V27ra s iP(0) B MU)j/ipC A C F (m b + m c f' 2 



N c m 3 J 2 (ml + m 2 — 2m c m b — q 2 ) 2 
16\/27ra s ^(0) 3^(0) j/^C A C Fy /m b + m c (m\ + Qm c ml + hm 2 c m b - q 2 m b + Ami ~ ^ m cq 2 ) 



a Lo _ j»v^"»fi»j5 e yi"j^ 

N c m Z J 2 (ml + m 2 — 2m c m b — q 2 ) 2 

(10) 

LO _ 32^^(0)^(0) .,/^C A C F (m b + m c fl 2 

a+ AT 3/2/ 2 , 2 o 2N2 ' ^ ' 

JS c m c [ml + — 2m c m b — q z y 

aL o _ 32y/27ra s ij(0) Bc ^(0)j/^C A C F y/m b + m c 
N c ^fm~ c (ml + m 2 c — 2m c m b — q 2 ) 2 

Here, the momentum transfer q = P — P', and the invariant mass q 2 — > 0, i.e. the finial 
charmonium owns the maximal momentum, denotes the maximum recoil point. 
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FIG. 2: The typical Feynman diagrams at one-loop level. 



III. THE NEXT-TO-LEADING ORDER CORRECTIONS 



In performing the next-to-leading order calculation, as schematically shown in Figure |2j 
we use the dimensional regularization scheme to regularize the UV and IR divergences, and 
the Coulomb divergence is regularized by the relative velocity v. In dimensional regulariza- 
tion, it is well-known that the 75 is difficult to deal with. In the literature, two approaches 
are mostly employed, that is the Naive scheme [l3] and the 't Hooft-Veltman scheme [l4]. In 
this calculation, we adopt the Naive scheme, of which the 75 anticommutates with each 7^ 
matrix in d-dimension space-time, £75, 7 M } = 0. In evaluating the quarkonium production 
and decays, it was argued by Ref. 15j that both schemes may lead to the same result, which 
is different from the case of pion decays to di-photon. The UV divergences exist merely in 
self-energy and triangle diagrams, which can be renormalized by the corresponding counter 
terms. The renormalization constants include Z2, Z3, Z m , and Z g , referring to quark field, 
gluon field, quark mass, and strong coupling constant a s , respectively. In our calculation the 
Z g is defined in the modified-minimal-subtraction (MS) scheme, while for the other three 
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FIG. 3: The ratios of NLO and LO form factors vs the square root of momentum transfer \/q 2 . 

pNLO/'2\ 

Here, Rj(<? ) = pzu^j with Fi standing for /+, /_, g, ao, a+, and a_. The renormalization-scale 
is fixed at = 3 GeV; = 4.76 GeV and m c = 1.54 GeV is adopted. 



fi 



the on-shell (OS) scheme is employed, which tells 

1 , 4tT// 2 



5Z° S 



5Z° S 



Q'., 



-a 



PCs 

An 



1 



7s + In —^r + - + 0(e 



3 

47T/X 2 



47T 



(A>-2CU)( 



1 



3 7B + 31n^- + 4 + C(e) 



+ 0(c) 



5Z 



MS 



' 2 An 



-7 S + ln47r + C(e) 



(13) 



Here, (3 = (11/3)Ca — (4/3)T/n/ is the one-loop coefficient of the QCD beta function; 
rif = 3 is the number of active quarks in our calculation; Ca = 3 and Tp — 1/2 attribute to 
the SU(3) group; \i is the renormalization scale. 

Because in our calculation the m c /mb contribution is kept, the complete expression turns 
to be too lengthy to be presented here. Whereas, the asymptotic form in small m c limit is 
given in the appendix, and the numerical results are presented. From the asymptotic form 
a noteworthy finding is that there exists an interesting relation among those form factors, 
i.e, 



,NLO 

l _0 

„LO 



,NLO 



9 



NLO 



9 



LO 



(14) 



which is consistent with the prediction of Ref.fl6| from large energy effective theory (LEET). 
In numerical calculation, the input heavy quark masses are 



m b 



4.76 GeV, m c 



1.54 GeV . 



(15) 



The one loop result for strong coupling constant, the 
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(11 — |n/)Log 



(16) 



V QCD 



is used. 



In Figure El the ratios of NLO and LO form factors versus the square root of momentum 
transfer y/q^ are schematically shown, while it should be noted that the pQCD approach is 
feasible only in the maximum recoil region. The figure shows that the NLO corrections to 
the B c to charmonia transition form factors are remarkable and sensitive to the momentum 
transfer q. The renormalization-scale dependence of the LO and NLO form factors are 
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FIG. 4: The renormalization-scale dependence of the LO and NLO form factors at the maximum 
recoil point, q 2 = 0. Here, Sj(/u) = j^ff ) with F{ standing for f + , /_, g, ao, a + , and a_. The 
solid line represents the situation of the LO renormalization-scale dependence and the dash line 
represents for the NLO situation. In the computation, m& = 4.76 GeV and m c = 1.54 GeV are 
adopted. 

presented in Figure H] at the maximum recoil point q 2 = 0. Generally speaking, the scale 
dependence in NLO is obviously depressed relative to the LO case. 

To show more explicitly the difference of LO and NLO results, we employ the function of 
SrXaO = 477 pLQ 1 } ) t° exhibit the dependence of renormalization-scale at the maximum 
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FIG. 5: The renormalization-scale dependence of the LO and NLO form factors at the maximum 
recoil point, q 2 = 0. Here, Sf^(/i) = 4na s Slof , with standing for / + , /_, g, ao, a+, and a_. 
The symbol "LO" denotes the LO renormalization-scale running of S^(^). In the computation, 
nib = 4.76 GeV and m c = 1.54 GeV are adopted. 




FIG. 6: The renormalization-scale dependence of the LO and NLO form factors at the maximum 
recoil point, q 2 = 0. In the computation, nib = 5 GeV and m c = 0.3 GeV are adopted and S^(/u) 
is performed to all orders of m c / nib- 
recoil q 2 = 0, as shown in Figures [5] tod where Fi(fi) stands for Fj i0 (/i) or F i LO (/i)-|-F i 7VLO (/x) 
corresponding to LO renormalization-scale dependence and NLO one respectively. To test 
the validity of m c / nib expansion, giving m c a nonphysical small mass we present in Figures 
|6]and [7] the leading order and full order results, respectively. From these figures we see that 
those form factors are quite sensitive to the magnitude of m c , and the higher order effects 
in m c /m b expansion are notable. 
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FIG. 7: The renormalization-scale dependence of the LO and NLO form factors at the maximum 
recoil point, q 2 = 0. The left one is for = 4.76 GeV and m c = 1.54 GeV; the right one for 
m;, = 5 GeV and m c = 0.3 GeV; and the S^(/j) is performed to the leading order of rric/mi,. 



grams, the FeynCalc 



In our calculation the package FeynArts [17j was used to generate the Feynman dia- 



18| was used to generated the amplitudes, and the Loop Tools [19] was 



employed to calculate the Passarino-Veltman integrals. 



IV. WAVE FUNCTION OVERLAP CONTRIBUTION 



The wave function overlap contribution for B c decays has been broadly discussed in the 



literature 



20 



331 ] . Since in the overlap contribution the QCD non-perturbative effects tend to 



be important, from pQCD factorization point of view it is beyond the scope of our interest in 
this work. However, here we still make a schematic comparison of the wave-function overlap 
contribution with pQCD contribution to the form factors for readers convenience. 

From the Table I, we notice that the results from wave function overlap are comparable to 
what from the pQCD calculation at the maximum recoil point, though in fact they attribute 
to different mechanisms in the study of B c decays. 
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TABLE I: B c 
overlap [33]. 



and J/V> form factors at maximum recoil q 2 = from pQCD and wave function 





| pB c r]c 1 


. R J lih i 
\ A Q \ 


, aB 7 /ib, 
| j^cj 1 iy | 


. B T/ib , 
\ A 2 \ 




DW[20] a 


0.420 


0.408 


0.416 


0.431 


0.591 


CNP[21] 


0.20 


0.26 


0.27 


0.28 


0.38 


KT[22J 


0.23 


0.21 


0.21 


0.23 


0.33 


KL0[23J 6 


0.66 


0.60 


0.63 


0.69 


1.03 


NW[24J 


0.5359 


0.532 


0.524 


0.509 


0.736 


IKS [25] 


0.76 


0.69 


0.68 


0.66 


0.96 


Kiselev[26j c 


0.66[0.7] 0.60[0.66] 0.63[0.66] 0.69[0.66] 1.03(0.94] 


EFG[27] 


0.47 


0.40 


0.50 


0.73 


0.49 


IKS2[28] 


0.61 


0.57 


0.56 


0.54 


0.83 


H1N V [29\ 


0.49 


0.45 


0.49 


0.56 


0.61 


HZ [30] 


0.87 


0.27 


0.75 


1.69 


1.69 


SDY[31J 


0.87 


0.27 


0.75 


1.69 


1.69 


DSV[32] 


0.58 


0.58 


0.63 


0.74 


0.91 


WSL[33] 


0.61 


0.53 


0.50 


0.44 


0.74 


pQCD(LO) a! 


0.97 


0.86 


0.90 


0.97 


0.32 


pQCD(NLO) d 


1.57 


1.34 


1.35 


1.39 


0.52 



°We quote the results with to = 0.6 GeV. 

b We quote the values where the Coulomb corrections are taken into account. 
c The results out (in) the brackets are evaluated in sum rules (potential model). 
d i[> Bc {0) = 0.36 GeV 3/2 , i/>j/V>( ) = °- 26 GeV 3/2 , a s = 0.2 



V. SUMMARY AND CONCLUSIONS 



In this work we have calculated the transition form factors of the B c meson to S-wave 
charmonia in the degree of NLO accuracy of pQCD. In our calculation, the B c meson and 
charmonia are treated as non-relativistic bound states of two heavy quarks. Hence, the 
long-distance effects for the form factors come only from the soft gluon exchange between 
heavy quarks, which can be explicitly factorized out and expressed as the wave functions at 



1 1 



the origin. The factorization scale is set to be at vm c in practical calculation, although the 
scale dependence does not appear at the one-loop level. Since in the small recoil region the 
end-point divergence spoils the QCD factorization, our result is valid only in the large recoil 
case. 

Calculation shows that the NLO QCD corrections to the B c to charmonia form factors are 
remarkable, especially for the /_ where the correction is as large as 80% or so. We find that 
the renormalization-scale dependence of the form factors are depressed, as it should be, when 
the next-to- leading order correction is taken into account, which means the uncertainties in 
the theoretical estimation are reduced. We find that B c to charmonia transition form factors 
are quite sensitive to the magnitude of m c , and the higher order effects in ra c /m& expansion 
are indispensable. 

Last, it should be mentioned that the relativistic corrections are also important for the 
B c to charmonia form factors, which are highly related to the relative velocity of heavy 
quarks within the bound states, i.e. ~ 0.38 and v\ c ~ 0.25 for B c and charmonium, 
respectively. 
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Appendix 

Following, various B c to charmonium transition form factors are given in leading power 

2 2 2 

of m c /mb and NLO pQCD. For the sake of compactness, we define s = m ^ 2 and 7 = ^ ^ . 

jNLO rNLO b 

It is worth emphasizing that our expressions for ~t L0 and —fut agree with what given in 
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reference 12J. 
fi (q 2 ) 



10n, 



7T 



61og(2))(s-l) + 3sbg( T ) 



9 6s + 3 

(18s 2 (2s - 1) log 2 (s) + 18(8 log(2)s 3 - 2 log(2)s 2 - 5 log(2)s + s 



72s 2 - 1 

+2 log(2)) log(s) + (2s - l)(268s + vr 2 (6s 2 - 3s - 6) + 170) - 9(2s 
-1) log( 7 )(log( 7 )s - (2 + 2 log(2))s + 41og(2)) + 18(2s - l)(4s 2 + s 
-2)Li 2 (l - 2s) - 18(4s 3 - 5s + 2)Li 2 (l - s) + 18(s(4s(s + 1) - 11) 
+4) log 2 (2) - 36(5(s - l)s + 1) log(2)) 



+ (5s - 2) log( 7 ) - 4(s + 1) log(2) + 12)(1 - 2s) 2 - 12(4s 2 + s 
-2)Li 2 (l - 2s)(l - 2s) 2 + 12(s(2s + 3) - l)Li 2 (l - s)(l - 2s) 2 
-(tt - 2vrs) 2 (s(4s - 19) + 4) + 3(-32 log 2 (2)s 4 - 4(69 + 21og(2)(-37 
+5 log(2)))s 3 + 8(18 + log(2)(-31 + 9 log(2)))s 2 + (61 + 28 log(2) 
-261og 2 (2))s + 121og(2) + 21og 2 (2) - 32) + (6s(8s(s(-41og(2)s 
+3 log(2) + 3) + 2 log(2) - 3) - 18 log(2) + 7) + 24 log(2)) log(s))} 



(17) 



/r°(Q) 

/f°(0) 



1 + — { 

4vr 3 



1 -(llC A -2n f ) log 
3 \m b m c J 



lOn, 



1 , ( m b 

3 V m. 



I. ',••*('** 

2 V m. 



- t !og( 2 ) log — 
3 \m r 



35 / m 6 
+ — log 



6 



+ 



21og(2) 
3 

21og 2 (2) 



, , N 7vr 2 103\ 

+C ,(-ilog=(^) + l log(2)log (M + l lof 
V 6 \m c J 3 \ m c/ 3 

41og(2) 5tt 2 73 

3 "36~ + ~9~ 



rn, 



— + 



log 2 (2) 



(18) 



13 



1 + ^ { 1 {UCa - 2n '> l0g( 2^ } ~ ^ + ^ bg(7) " 6 l0g(2) + " 2) 

+ 36(s -l) A (2s - l) (18(g " 1)S(2S " 1} l0g2(s) + 18 ^ 2s ( 4s " 5 ) + X ) M 2 ) s 
+s + log(2) + 1) log(s) + (s - l)(2s - 1)(tt 2 (6s - 3) + 268) 

+9(s - l)(2s - 1)(- log( 7 ) + 21og(2) + 2) log( 7 ) + 18(8s 3 - 10s 2 + s 

+l)Li 2 (l - 2s) - 18(s - l)(2s - l)(2s + l)Li 2 (l - s) + 18(4s 3 - 7s 

+3)log 2 (2)-36(s-l)(5s-l)log(2)) 

(12(1 - 2s) 2 log 2 (s)(s - l) 3 + (tt 2 (s - l)(4s - 19)(1 - 2s) 2 



12(l-2s) 2 (s- 1) 

+3(32 log 2 (2)s 4 + 4(69 + 2(-37 + log(2)) log(2))s 3 + (-508 - 8 log(2)(-74 
+13 log(2)))s 2 + (307 - 364 log(2) + 82 log 2 (2))s + 2(34 - 9 log(2)) log(2) 
-61))(s - 1) + 6(s(s(-24s 2 + 84s + 2(2s - l)(2s(4s - 9) + 11) log(2) - 127) 
-4 log(2) + 73) + 2 log(2) - 13) log(s) + 3(2s 2 - 3s + l) 2 (-5 log( 7 ) + 4 log(2) 
-23) log( 7 ) + 12(4s + l)(2s 2 - 3s + l) 2 Li 2 (l - 2s) 

-12(2s + 3)(2s 2 - 3s + l) 2 Li 2 (l - s)) (19) 

1 + S ^_ 2ny)l0g (^)_^4 l0g g)_ 2l0g(2) 4 

+ lllog(2) + 5^ _ \ 
2 3 J 

+C A (-i-vj^+l log(2) log (Si) + I , og (H) + log2 < 2 » 



4 \m c J 2 \m c J 2 \m c J 2 

7T 2 76 \ 

-51og(2)- y + ^j (20) 



14 



" 1 + 4^ { 3 (11 ^- 2n/)1 ° g( 2^ ) -^ 
C 

1 : (9s(2s-l)log 2 (s) + 18(2slog(2)(2s-l) + l)log(s) 



36s- 18 

+3vr 2 (s + 2)(2s - 1) - 2s(-18 log 2 (2)s + 9 log 2 (2) + 45 log(2) + 134) 
+9(2s - l)(log( 7 ) - 3) log( 7 ) + 18s(2s - l)(2Li 2 (l - 2s) - Li 2 (l - s)) 
+63 log(2) + 134) 

6(1 _ 2 ^ 2(g _ 1} (6( g 2 - 1) log 2 (s)(l - 2s) 2 + 24(s - l)sLi 2 (l - 2s)(l - 2s) 2 
+3(2s(s(4s(4 log(2)s - 8 log(2) + 3) + 20 log(2) - 17) - 4 log(2) + 7) - 1) log(s) 
+ (s - 1)(6 log( 7 )(log( 7 ) - 6 log(2) + 5)(1 - 2s) 2 + 6(2s - 9) log 2 (2)(l - 2s) 2 
+ (2s - l)(-204s + 2tt 2 (2s 2 + s - 1) + 105) + 6(s(68s - 67) + 16) log(2)) 
-12(2s 2 -3s + l) 2 Li 2 (l-s))} (21) 



|i { i(llC,-2, V )log(^l) 

47T 3 \m b m c J 



g NLO (0) a Sf l /iiri o Nl ( 2/i 2 ^ 10n/ 



+C F ( log 2 ( ^ ) - 10 log(2) log ( ^ ) + 5 log ( ^ ) + 9 log 2 (2) 



m c J \m r I \m r 

+71og(2) + y-15 

31og(2) ^ + 67A 



2 3 9 



i$ LO (q 2 ) a N + LO {q 2 ) g NLO (q 2 ) 



a%°(q 2 ) a^iq 2 ) g LO (q 2 ) 

a$ LO (0) _ < LO (0) _ g NLO (0) 
a^°(0) ~ a£°(0) ~ g LO {0) 



(23) 



(24) 
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a 



NL °tf) . 1 , <*' f l,„n o M , 10n / 



1 + 7zi o(^C A -2n f ) log 



a LO (g 2 ) 4tt 1 3 v j; 2^2/ 9 

-(18(2 log(2)s 2 - 3 log(2)s + s + log(2)) log(s) 



9(s- l)(2s- 1) 

+ (s - l)(2(-85 + 36(-l + log(2)) log(2))s + vr 2 (6s - 3) + 18(2s - 1) log(2) log( 7 ) 

-181og(2)(-l + 21og(2)) + 85) + 18(s - l)(2s - l)Li 2 (l - 2s) 

-18(s-l)(2s-l)Li 2 (l-s)) 
C F 



(log 2 (s)(2s 2 - 3s + I) 2 + 4Li 2 (l - 2s)(2s 2 - 3s + iy 



(2s 2 -3s + l) 5 

-2Li 2 (l - s)(2s 2 - 3s + l) 2 + (s(s(4(3 + 41og(2))s 2 -8(5 + 6 log(2))s 

+52 log(2) + 27) -4(1 + 6 log(2))) + 4 log(2)) log(s) + |(12(-13 

+2 log(2) + log 2 (2))s 4 - 12(-43 + 5 log(2) + 3 log 2 (2))s 3 + (-597 

+54 log(2) + 26 log(2) 2 )s 2 - 3(-95 + 6 log 2 (2) + 8 log(2))s 

-2tt 2 (2s 2 - 3s + l) 2 - 3(2s 2 - 3s + l) 2 log( 7 )(log( 7 ) + 2 log(2) - 6) 

+31og 2 (2) + 61og(2)-48)) (25) 



«. 1 + ^ (11Cj _ 2n/)log (^)_^ 

47r 3 \mbm c J 9 



+C F (■ log 2 [^-J +21og(2)log • 6 log ) bg 2 (2) 

-61og(2)-7r 2 -^ 
+CU (-2 log(2) log + 6 log(2) - y + f ) (26) 
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